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We have studied effects of interspecies attraction in a Fermi-Bose mixture over a large regime of 
particle numbers in the / ^^Rb system. We report on the observation of a mean field driven 
collapse at critical particle numbers of 1.2 • 10^ ^^Rb atoms in the condensate and 7.5 • 10^ '^^K 
atoms consistent with mean field theory for a scattering length of a^s = —281(15) ao [S. Inouye et 
al, Phys. Rev. Lett. 93, 183201 (2004)]. For large overcritical particle numbers, we see evidence 
for revivals of the collapse. Part of our detailed study of the decay dynamics and mechanisms is a 
measurement of the (^'^Rb - ^'^Rb - ^°K) three-body loss coefficient K3 = (2.8 ± 1.1) • 10"^^cmVs, 
which is an important input parameter for dynamical studies of the system. 
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The recent realization of the BCS-BEC-crossover 
in strongly interacting dilute fermionic gases has allowed 
intriguing insight into the regime connecting Bose and 
Fermi superfluidity. Quantum degenerate Fermi-Bose 
mixtures are expected to offer an alternative and comple- 
mentary approach to Fermi superfluidity, where the in- 
teraction between fermions is mediated by the bosons Q , 
analogous to the role of phonons in solid state supercon- 
ductors. The realization of boson-mediated Cooper pair- 
ing in optical lattices with tunable interactions through 
hetero-nuclear Feshbach resonances would open up fas- 
cinating perspectives. Furthermore, it has been pointed 
out that these mixtures exhibit a wealth of novel quan- 
tum phases 0. 

Exploring these phenomena relies on a detailed un- 
derstanding of interactions in ultracold Fermi-Bose mix- 
tures which are characterized by two scattering param- 
eters: the s-wave background scattering length for colli- 
sions between two bosons asB and the Fermi-Bose inter- 
species scattering length apB- For repulsive interspecies 
interaction, the constituents may phase separate, while 
an attractive interaction creates an additional mean field 
confinement and eventually renders the mixture unsta- 
ble 0,1^, similar to the mean field collapse found in BECs 
with attractive interaction 6]. 

The prediction and first experimental confirmation of 
a large and negative interspecies scattering length 0, 
shifted the ^^K / ^^Rb system into the focus of scien- 
tific interest. In addition to the very precisely known 
dBB = 98.98(04) ao 0, the value oiapB is of paramount 
importance for this system, setting maximum achievable 
particle numbers but also determining correlations and 
excitations. Published critical particle numbers for the 
occurrence of a mean field collapse {Nk ^2-10^ and 
Nuiy ^ 10^ ^3) appeared to impose severe constraints 
on maximally achievable particle numbers, but at the 
same time seemed to indicate an excitingly large value 
oi apB = -395(15) ao There is however an on- 

going controversy in the scientific community caused by 
recent reports on the observation of stable mixtures at 



higher particle numbers ^|, and measurements of 
apB = — 250(30) ao using thermal relaxation and 
aFB = — 281(15) ao |l3 using Feshbach spectroscopy 
methods. A detailed exploration of the regime of mean 
field instability is not only interesting in itself due to the 
its strong dependence on the interaction parameters and 
due to the dynamical behavior when the instability oc- 
curs. More generally, it is a fundamental prerequisite for 
a full understanding of ^^K / ^^Rb-mixtures and the con- 
trolled achievement of fascinating physics in the regime 
of large particle numbers and strong interactions, e.g. the 
realization of mixed bright solitons Q and the observa- 
tion of Fermi-Bose correlations up to novel BCS phases 
in optical lattices. 

In this letter, we report on reaching the to our knowl- 
edge so far highest particle numbers in the ^^K / ^^Rb 
system only limited by the onset of instability, which we 
find consistent with apB ^ — 281ao 0. Right at the 
phase transition point of a large bosonic thermal cloud 
(n ~ 4 • 10^^ cm~^), we observe a strong dynamical mod- 
ification of the in-trap fermionic distribution due to the 
mean- field interaction with the bosons. At even higher 
densities, we observe enhanced localized loss processes in 
the overlap region of the BEC and the Fermi gas, which 
we ascribe to a mean field collapse based on lifetime and 
dynamical arguments as well as comparison to theory. 
After the collapse, we observe a peaked density distribu- 
tion of the Fermi cloud subjected to the strong localized 
mean field potential of the condensate. In addition, we 
present a novel measurement of the three-body inelastic 
K-Rb decay constant differing by an order of magnitude 
from the value given in [lo| . 

Our setup is based on a two-species 2D/3D-MOT ge- 
ometry 15], where the ^^K vapor is produced using self- 
made enriched ^^K dispensers In the 3D-M0T, we 
achieve 1-10^^ ^^Rb and 2-10^ fermionic ^^K atoms within 
10 seconds. The pre-cooled mixture is loaded into a loffe- 
Pritchard type magnetic trap (cloverleaf / 4D-hybrid, 
uf^ = 11.3 Hz and u^J^^ = 257 Hz, i.e. uf^ = 16.6 Hz 



and 



378 Hz) for RF-induced sympathetic cooling 
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of ^^K in the |F = 9/2, mF = 9/2 > state with ^^Rb in 
the |F = 2, mp = 2 > state. If all of the ^^Rb atoms are 
removed, we currently obtain Fermi gases of up to 3 • 10^ 
^^K atoms at T/Tf = 0.2 or 9 • 10^ at T/Tf = 0.1. 

When the ^''Rb atoms are not completely evaporated, 
various regimes of mixtures are accessible, ranging from 
dense thermal ^''Rb clouds of 10^ ^''Rb atoms right at 
the phase transition point interacting with a moderately 
degenerate Fermi gas of 2-10^ "^^K atoms to deeply degen- 
erate mixtures with almost pure condensates. We achieve 
> 1 • 10^ atoms in the condensate coexisting with 7.5 • 10^ 
^^K atoms, limited by the onset of the collapse. 

In a first set of experiments, shown in Fig. ^ we ana- 
lyze the behavior of the mixture in these regimes during 
the evaporation ramp, similar to the procedure in In 
our case the dynamics becomes clearly visible in the ax- 
ial ^^K density distribution, which due to the relatively 
short time of flight (TOF) should closely reflect the in 
trap distribution. The fast ramp speed (IMHz/s) in 
this experiment means that we are potentially dealing 
with strongly out of equilibrium samples. Already at the 
phase transition point (a) of the bosonic cloud, we ob- 
serve strong distortions of the axial "^^K density profile. 
The profile looks like a chopped off Fermi profile with 
a peak in the center of the flat top which we ascribe to 
the interaction with the bosonic component (peak den- 
sity ^ 4 • 10^^ cm-3). For a large BEC (b) the "^^K profile 
exhibits a pronounced hole in the trap center which we 
ascribe to a strong localized loss process due to the inter- 
action with the BEC. These losses are too fast for trans- 
port in the slow axial direction to be able to continuously 
maintain the undisturbed density profile by refilling the 
center of the trap from the outer regions. As we shall see, 
this fast loss (e. g. from Fig. [2l < 10 ms) is a signature 
of the mean field collapse of the mixture. After the col- 
lapse, we observe that ^^K distributions sharply peaked 
in the center remain stable for relatively long timescales 
exceeding 100 ms (c). A Fermi-Dirac fit to the data is 
shown as a dotted line for comparison. We attribute this 
peaked distribution to the Fermi-Bose attraction creating 
an additional trapping potential for the fermions ( "mean 
field dimple" in the magnetic trap potential). 

The strong depletion of the Fermi cloud in the center 
is accompanied by rapid particle loss of approximately 
two thirds of the fermionic particle number. We ascribe 
this loss to the mean field collapse of the mixture. Be- 
yond critical particle numbers, the mean field potential 
is no longer balanced by the repulsive interaction in the 
^''Rb BEC and the outward bound Fermi pressure, so 
that the part of the mixture overlapping with the BEC 
contracts rapidly to such large densities that enormous 3- 
body losses reduce the overall particle number in this re- 
gion to an under critical value. It is therefore natural that 
we observe the collapse in the vicinity of strong contract- 
ing mean field effects of the Bose gas on the fermionic 
distribution (see Fig. c; cf. ^3)- Although the en- 




axial coordinate [arb. un.] 



FIG. 1: Typical evolution of an overcritical mixture after an 
evaporation ramp (rate -1 MHz/s) stopped and held fixed for 
15 ms at 80 kHz (a), 50 kHz (b) and 20 kHz (c) above the ^'^Rb 
trap bottom of 490 kHz. Left-hand side: 3D representation 
of absorption images with false-color coding of the optical 
density. ^^Rb and images are taken in the same run, 
although at different TOF: 20ms (^'^Rb) and 3 - 5ms (^°K). 
Right-hand side: corresponding axial line profiles integrated 
along the vertical direction. 



semble is out of equilibrium and shows very complex dy- 
namics after this rapid loss , one can intuitively imag- 
ine that the Fermi distribution depleted in the center of 
the cloud is refilled from the outside parts of the sample 
on a timescale related to the axial trap frequency, possi- 
bly leading to repeated local collapses, until the mixture 
will become undercritical and reestablish an equilibrium 
situation. In order to observe this phenomenon we have 
prepared a mixture where the bosonic part (~ 10^ atoms) 
has only started condensing and observed the evolution 
of the ensemble at constant evaporation frequency. As 
the condensate grows due to the mild remaining evapo- 
rative cooling, it reaches the critical particle number of 
Nb = 1.2 • 10^. Due to the near-equilibrium situation 
before, the collapse now only leads to a relatively small 
loss in total particle number but is still clearly visible in 
the ^^K atom number integrated over the central part of 
the TOF image as shown in Fig. [2j After the first col- 
lapse this number remains constant for some time and 
then drops abruptly again 26]. Such "revivals" of the 
collapse have been predicted in a recent numerical anal- 
ysis of the collapse dynamics (^^, see inset in Fig. [SJ, 
although in a spherically symmetric configuration with 
i/Rb = 100 Hz. 

Analysis of similar decay series in various particle num- 
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FIG. 2: Decay of the fermionic component in a slightly over- 
critical mixture. The first sudden drop is the initial collapse; 
the second drop is one "revival" of the collapse. The line is to 
guide the eye. For comparison, the inset shows results from 
dynamical modelling of the collapse. 




FIG. 3: Stability diagram for the / ^'^Rb mixture. The 
uncertainty in A^Rb/A^K is assumed to be 20%/30%. The error 
on Nf is not critical and is given exemplarily in (b). The 
solid black line is based on the theory of ref. ^3 clfb = 
-281(15)ao 111. (a) is the critical particle number reported 
in ref. [ifl (here the trap had a different aspect ratio, but a 
similar mean trapping frequency i>K = 134 Hz, as compared 
to our experiment with i>K = 133 Hz) 



ber regimes enables us to extract a value for the criti- 
cal particle numbers for the onset of the collapse. Our 
findings are summarized in Fig. |31 Situations where 
the decay is compatible with 3-body decay of an undis- 
turbed density distribution are identified as stable. If 
the loss is incompatible with 3-body decay, the situations 
are identified as unstable. Note that the time scales for 
the two situations are clearly different, varying from a 
few 100 ms for the stable cases to about 20 ms for the 



unstable cases. The solid line in Fig. |31 is the theo- 
retical stability limit for our trapping conditions, based 
on 19] and the following approximation to the fermionic 
chemical potential: /ii? = /ii?,o — {sfb/ Sbb) ' Mb,o where 
fiF,o and /iB,o are the Fermion / Boson non-interacting 
chemical potentials an d qj ^b and gsB are the standard 
coupling constants |20, l2l| . The curve is plotted for the 
value aFB = — 281ao; a confidence interval for the stabil- 
ity border based on the uncertainty of 15ao on apB 
is indicated in grey. Our findings for the critical particle 
numbers are clearly consistent with mean-field theory of 
the collapse based on this value for apB ^27]. 

Three-body decay is the underlying loss mechanism for 
the collapse of the mixture. The value of the associated 
rate coefficient is therefore an important input parame- 
ter for detailed understanding of the dynamical behavior 
of the collapse ("revivals") 18]. In order to minimize 
experimental uncertainties, we extract the value for the 
3-body rate coefficient from the decay of non- degenerate 
thermal ensembles in an equilibrium situation over a large 
range of densities. Since collisions involving two identical 
fermions are strongly suppressed due to the Pauli exclu- 
sion principle US, 13, we consider only inelastic collisions 
between two bosons and one fermion. For pure 3-body 
loss, the decay process is then characterized by the fol- 
lowing rate equation: 
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where ub and np are the respective densities, is the 
rate coefficient for three-body loss and r the background 
collisional lifetime. Integrating this equation over time 
yields 
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Based on the measured background scattering rate, the 
left-hand side expression can be evaluated at any given 
time T for a decay series. The right-hand side inner 
integral is evaluated using an equilibrium Fermi-Dirac 
profile for ^^K and a Bose-enhanced thermal profile for 
^^Rb at different time values ti . The density distributions 
are based on the measured temperatures and total atom 
numbers at t^. The outer integral is evaluated by using 
a piecewise linear interpolation between the discrete val- 
ues of the inner integral. Plotting the right-hand side as 
a function of the double integral, we obtain i^s as the 
slope of the linear plot (see inset (a) in figure^. Note 
that our measurement does not depend on the potas- 
sium atom number calibration. As far as the ^^Rb atom 
number calibration is concerned, we reproduce the mea- 
surement of ^^Rb 3-body decay from ref. l22l . We find 
a value of = (2.8 ± 1.1) • 10"^^ cm^/s [ij nearly an 
order of magnitude lower than the value found in [lOj . 

The additional analysis of decay in degenerate samples 
at higher densities reproduces this result as long as the 
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FIG. 4: Extraction of 3-body inelastic loss rate, (a) decay 
analysis of a thermal sample, where J\f{t) = In j^^^^^ + ^ and 

JC{t) = N-^ . jldt' j dPrnl{r,t') ■ nF{r,t') (b) ii^s values ex- 
tracted in the moderately degenerate and stable regime with 
peak BEG densities of up to 4 • 10^^ cm~^ are compatible with 
the thermal average (solid line), (c) For low evaporation 
end frequencies and BEG densities of ^ 5 — 7 • lO^^cm"^, 
samples are unstable with respect to the collapse. 



mixture is stable (03), but shows significant deviations 
for unstable mixtures (^). In the regime of instability, 
our analysis of the loss process is no longer valid. Due to 
the strongly depleted ^^K density distribution, our analy- 
sis, relying on undisturbed density distributions, overesti- 
mates the density overlap integrals which should lead to a 
too low value of as extracted from the decay. Instead, 
as can be seen from Fig. increases below evapo- 

ration end frequencies of 60 kHz (high initial densities), 
proving that the underlying very fast loss mechanism is 
due to the collapse. 

In conclusion, we have presented measurements cover- 
ing the so far highest particle numbers available in quan- 
tum degenerate "^^K-^^Rb Fermi-Bose mixtures leading 
from stability to collapse. Our measurements of critical 
particle numbers are consistent with a mean field model 
of the collapse based on the value of the scattering length 
apB = — 281(15) ao obtained from Feshbach resonance 
data This settles the clear discrepancy between the 
theoretical description of the collapse based on the above 
interspecies scattering parameter and the critical parti- 
cle numbers for the mean field collapse reported in an 
earlier experiment (see Fig. and indicates that 
these findings might be influenced by some additional 
interesting effect. Part of our detailed study of the mix- 
ture is a quantitative analysis of three-body loss leading 
to an order of magnitude lower value than in (i^. Our 
findings indicate that a wide range of filling factors in 
large- volume optical lattices ranging in principle up to 2 
are possible for the "^^K-^^Rb system. Progress towards 
loading the mixture into a 3D optical lattice is currently 



under way in several labs. Finally, reaching the regime 
of mean field collapse also lays the ground for the re- 
alization of bright solitonlike structures in a quasi ID 
geometry 
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